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Abstract 

The method of In bilayer pre-deposition and penetrated nitridation had been proposed, which had been proven to 
have many advantages theoretically. To study the growth behavior of this method experimentally, various pulse 
times of trimethylindium supply were used to get the optimal indium bilayer controlling by metalorganic vapour 
phase epitaxy. The results revealed that the InN film quality became better as the thickness of the top indium 
atomic layers was close to bilayer. A following tuning of nitridation process enhanced the quality of InN film further, 
which means that a moderate, stable, and slow nitridation process by NH3 flow also plays the key role in growing 
better-quality InN film. Meanwhile, the biaxial strain of InN film was gradually relaxing when the flatness was 
increasingly improved. 
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Background 

The unique properties of InN are currently attracting 
much interest in the research community [1,2]. Because 
of its lowest effective mass and the highest electron drift 
velocity among all Ill-nitride semiconductors [3], InN is 
promising for high-speed and high-frequency electronic 
devices. And recently, the band gap of InN, which is 
considered as 1.9 eV, is renewed to approximately 0.7 eV 
[4-6], covering a broad range of wavelength from near 
infrared at approximately 1.5 um to ultraviolet at ap- 
proximately 200 nm based on its direct band gap alloy- 
ing with GaN and A1N [7-9]. However, the achievement 
of high crystalline quality InN film is still a great chal- 
lenge, which attributed to the lack of lattice-matched 
substrate, the low dissociation temperature (approxi- 
mately 600°C), the low pyrolysis efficiency of NH 3 at 500 
to 600°C, the pre-reaction of the precursors before arriv- 
ing at the substrate surface, and also the InN clustering 
effect [10-14]. 

In order to achieve high-quality InN film, effort has 
been made by researchers with different methods such 
as optimizing growth temperature, controlling V/III ra- 
tio, introducing buffer layer, or employing pulsed atomic 
layer epitaxy technique [15,16]. However, the crystalline 
quality of InN film is still far below a satisfactory level 
due to the existence of huge quantity of defects [16]. To 
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elucidate the original difficulty in In film deposition, the 
formation kinetics of InN with N and In atoms on the 
In polar GaN surface has been systematically studied by 
first-principles calculations [17], it was found that the 
pre-deposition of In bilayer on the surface could im- 
prove the In migration on the surface and the smooth- 
ness of In film. 

In this work, the epitaxy method of In bilayer control- 
ling and penetrated nitridation was employed for the 
InN film growth on GaN template. In order to deter- 
mine critical trimethylindium (TMI) flow required for 
forming In bilayer, the pulse time of TMI supply was op- 
timized. The results revealed that the film quality be- 
came better as the thickness of the top indium atomic 
layers was close to bilayer. Based on the In bilayer de- 
position, a moderate, stable, and slow nitridation process 
by NH3 flow also played the key role in growing better- 
quality InN film. X-ray diffraction (XRD) measurements 
confirmed the gradual relaxation of biaxial strain in InN 
epilayers during increment of the smoothness. 



Methods 

Growth of samples 

InN films were grown on a 3-um-thick GaN template 
with(0001) sapphire substrate by using metalorganic 
chemical vapor deposition (MOCVD) system with a 
Thomas Swan closely coupled showerhead (CCS) re- 
actor. The trimethylgallium (TMG), trimethylindium 
(TMI), and ammonia (NH 3 ) were used as the precursors 
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for Ga, In, and N, respectively, and H 2 and N 2 were used 
as the carrier gasses. Prior to the GaN/AlGaN superlat- 
tice growth, thermal cleaning of the (0001) -oriented sap- 
phire substrate was carried out under hydrogen ambient 
at 1,050°C for 10 min to remove native oxide from the 
surface. Then, an approximately 30-nm low-temperature 
GaN buffer layer (approximately 570°C) was grown 
followed by a approximately 3-um high-quality GaN un- 
derlaying layer (approximately 1,090°C). During the stage 
of InN growth, the pressure was set to 450 torr at 550°C 
[18]. In order to accurately control the deposition of 
indium atomic multilayers and the following nitrida- 
tion process, the pulse growth method was employed 
through switching and adjusting the pulsed supply 
time of TMI and ammonia flows, as shown in Figure 1. 
For samples A, B, C, and D, a constant TMI flow of 
2.0 x 10~ 5 mol/min was used whereas a series of dur- 
ation time of the pulsed TMI flow, 16, 8, 4, and 3 s, 
was applied, respectively. Then, they were followed by 
a 33-s pulse of NH 3 flow for the nitridation process. 
The mole flow of ammonia was set to be 0.5 mol/min. 
For the purpose of comparison, the total moles of 
indium supply in different samples was preserved by 
varying the total pulse periods of 180, 360, 720, and 
960 for samples A to D, respectively. In order to avoid 
the possibility of over-nitridation, a lower ammonia 
flow of 0.25 and 0.125 mol/min was set for samples E 
and F, respectively, and other growth parameters were 
as the same as those set for sample C. 

Characterization 

The thickness of film was measured by in situ growth 
monitoring curves (Panalytical Xpert PRO X, Panalytical, 
Almelo, The Netherlands). The surface morphology and 
smoothness of the as-grown samples were characterized 



by atomic force microscopy (AFM, PicoSPM and PSI 
XE-100, Molecular Imaging, Ann Arbor, MI, USA) and 
scanning electron microscopy (SEM, LEO 1530, LEO 
Elektronenmikroskopie GmbH, Oberkochen, Germany) 
equipped with an energy-dispersive X-ray spectrometer 
(EDX). The structural quality and the In composition of 
InN films were evaluated by X-ray diffraction (XRD) in a X' 
Pert PRO system. 



Results and discussion 

In the ideal indium bilayer construction process, we 
need to deposit one indium monolayer in each pulse, 
thus this new indium monolayer would construct an in- 
dium bilayer with the top indium monolayer which we 
had deposited in last pulse [17]. Figure 2 shows the in 
situ growth monitoring interferometer curves of samples 
A to D. One can observe that in the InN growth stage, 
the vibration of all four sample's curves have experi- 
enced nearly equal phase shift. According to this phase 
shift, we can easily calculate the average InN film's thick- 
ness of them, which is about 170 nm. Also this result has 
been confirmed by practical measurement in the SEM 
cross-sectional observation (see SEM cross-sectional pho- 
tos in Figure 3). Thus, the InN deposition thicknesses 
per period of samples A, B, C, and D are about 9.5, 4.7, 
2.4, and 1.8 A, respectively. According to the (0001) 
lattice constant c of wurtzite InN, the thickness of one 
In-N monolayer (c/2) is about 2.8 A. Comparing to this 
value, sample C's growth thickness of each pulse is the 
closest one. 

Figure 3 shows the SEM images of surface morphology 
and cross sections of samples A to D. From the top view 
of sample A (Al), one can see some obvious dark holes 
on the surface, indicating the formation of vacancies 
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Figure 1 The pulse timing's schematic diagrams of samples A to D. 
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Figure 2 The in situ growth monitor interferometery curves of samples A to D. 



due to In accumulation in droplets. The formation of 
holes and droplets easily leads to a pretty rough surface 
(rms = 33, from AFM scanning result), as shown in 
Figure 3A2. As we know, the melting temperature of 
metal In is only about 157°C Thus, under the growth 
temperature of InN (550°C), the pulsed deposition of 
In for a long duration time may form a thick liquid In 
layer on the surface. By the effect of surface tension, In 
droplets in large size would come into being quickly. This 
is the main reason governing the surface roughness. In 
order to reduce the roughness, the pulse time of TMI 
is reduced to 8 s for sample B. The obtained InN film 
shows better flatness (rms = 20) and dark holes have 
been well removed (Figure 3B2). According to the the- 
oretical simulation of the kinetics of InN formation 
[17], if the thickness of indium film is larger than two 
atomic layers, the nitridation of this In film could not 
well form a InN epilayers in correct stoichiometric ratio 
(1:1) and the excessive In will lead to roughness. Thus, the 
TMI pulse time was further decreased down to 4 s. As 
shown in Figure 3C1, the islands of sample C begin to 
show regular shape relatively and the surface becomes 
more flat (rms = 14). Meanwhile, it can be observed that 
there are some islands in larger size, as indicated by arrow. 
The number of these types of large islands further in- 
creases in sample D (Figure 3D1), in which the TMI 
pulse time was set to 3 s. This trend of quality deterior- 
ation implies that the indium film deposited during the 
TMI period turns to be less than one atomic layer and 
fail to construct indium bilayer. This insufficient cover- 
age of indium layer could not provide the advantage of 



nitridation of indium bi-layer structure. On the contrary, 
over-nitridation under N-rich condition leads to the de- 
terioration of the InN film quality of sample D. Therefore, 
it could be determined that 4-s pulsed supply of TMI in 
sample C is the optimal setting. 

To investigate the optical property of these samples, 
absorption spectra were recorded to determine the band 
gap of InN film and the results are shown in Figure 4. 
Although all four samples' absorption curves show lim- 
ited differences due to the small thickness or relatively 
low crystalline quality of the InN film, their differences 
of slope's changes still can be identified. The absorption 
spectra of sample C and D have a clear slope threshold 
near the absorption edge. While, for samples A and B 
such slope threshold is absent and, beyond 1,100 nm, 
absorptions related with defect or impurity bands ap- 
pear. This indicates that sample C has the best film qual- 
ity due to the optimized pulsed growth with TMI supply. 
In principle, InN is a direct band semiconductor so that 
the relationship between its energy band gap and its ab- 
sorption coefficient could follow the formula below: 

(ahvf = k(hv-E g ), (1) 

where the a is the absorption coefficient and the E g is 
the band gap. Thus, the E g of our samples could be esti- 
mated through the intersectional point of absorption 
edge's tangent and horizontal axis. It is found that the E g 
of sample C and D is about 1.22 and 1.19 eV, respect- 
ively. Due to the unclear slope thresholds in samples A 
and B, the E g is difficult to determine precisely. The 
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Figure 3 SEM images of samples A to D. (A1, B1, CI D1) The top view and (A2, B2, C2, D2) the side view images of samples 
A to D, respectively. 



range of reasonable E g for samples A and B would be be- 
tween 0.7 to 0.9 eV, which is lower than those of sample 
C and D. This may stem from the absorption associated 
with the deep level defects and impurities rather than 
the band edge absorption of InN, the E g of which should 
be about 1.25 eV [19]. 



Considering the negative influence by the excessive 
NH 3 supply, we tried to improve the nitridation process 
by optimizing the ammonia flow. In principle, the in- 
dium bilayer will experience a nitridation process with 
the penetration of N atoms into between the bilayer 
[17]. This process would finally form a uniform wurtzite 
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Figure 4 The absorption spectra of samples A to D. 



InN structure on the surface. For the case of excessive 
NH3 flow, the top layer in high N concentration on the 
surface easily forms a steep concentration gradient be- 
tween surface and sub-surface layers where the N atoms 
will diffuse to. According to Fick' first law, 



-D 



dC 



(2) 



where the / is the total diffusion flux and the D is the 
diffusion factor. The steeper the concentration gradient 
(~ i§) wou ld lead to the higher the total diffusion flux / 
[20]. Thus, N atoms could not uniquely arrive at the 
preferable top site via the one-atom-on-one-site mode. 
Instead, they would diffuse to various positions and 
some would even crowd in some energy minima. Mean- 
while, ultra-high N concentration on surface could even 
make some N atoms hang over the top indium atomic 
layer, and, in this case, the indium pre-deposition of next 
pulse would fail to construct indium bilayer in some re- 
gions. As a result, the uniformity and smoothness of the 
InN film is deteriorated. Based on this analysis, the NH3 
flow should be optimized by reducing the mass flow, 
which is set to 0.25 mol/min for sample E and 
0.125 mol/min for sample F. Figure 5 shows the SEM 
images of these two samples. One can see that the 
smoothness of sample E has been slightly improved and 
is better than that of sample C. This indicates that the 
lower ammonia flow could improve the uniform diffu- 
sion of N atoms. Further reduction of NH3 flow in sam- 
ple F finally leads to a large improvement of InN quality 
and surface smoothness, as shown in the cross-sectional 
image of Figure 5F2. The corresponding AFM scanning 
also confirms this enhancement of surface smoothness 
(rms = 7). After the deposition of indium bilayer, a 
moderate, stable, and slow nitridation process with 





Figure 5 SEM images of sample E and F. (El, F1) The top view and (E2, F2) the side view images of samples E and F, respectively. 
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appropriate ammonia flow is crucial for the formation of 
better-quality InN film. 

In order to study the residual strain of as-grown InN 
films, XRD characterizations with G)-26 scans were taken 
and the results are shown in Figure 6. Typical symmet- 
rical (002) diffraction peaks of wurtzite InN and wurtz- 
ite GaN could be clearly identified, at about 15.8° and 
17.4° [21]. Besides, another weak peak was observed at 
about 16.65°; this peak has been identified as (101) dif- 
fraction peak of wurtzite InN by consulting related 
database and reference. In order to separate the mixing 
of these two peaks, a multi-peak fitting in this region 
was made and peak positions of each could be deter- 
mined. Base on angle of the (002) peaks (green dash 
lines) and angle of the (101) peaks (red dash lines), the 
corresponding lattice constant along c axis and a axis 
could be obtained simply by using the formula below: 




16.0 16.5 17.0 17.5 

Theta(°) 

Figure 6 The XRD diffraction spectra of samples A, B, C, E, 
and F. 



v. J 



where the d^ki is the interplanar distance between (hid) 
planes. The complete crystalline data is summarized in 
Table 1. One can see that the lattice constant a is in- 
creasing from samples A to F, and the a value of sample 
F (3.63 A) is very close to the equilibrium value of 
wurtzite InN (3.627 A) obtained by first principle calcu- 
lations, indicating the gradual reduction of residual bi- 
axial strains through growth optimization. Whereas, the 
(002) peak (correspond to lattice constant c) is right 
shifting correspondingly due to the expansion distortion 
by the elastic strain on the a axis. Meanwhile, it can be 
seen that the (002) peak is getting dominant, which 
means a preferential (002) crystal orientation in sample 
F. All these evidences imply that the biaxial strain has 
been well relaxed, and the crystal orientation has be- 
come better in sample F. 



Conclusions 

Through using various pulse times of TMI supply, 
we achieved optimal indium bilayer control by meta- 
lorganic vapour phase epitaxy. When the top indium 
multilayer was getting close to bilayer, InN film 
quality had been gradually improved due to high 
surface migration and good structure consistency of 
indium bilayer forming. The absorption spectra also 
confirmed that the InN film which was grown via opti- 
mal indium pre-deposited controlling had the fewest 
defects and impurities. Furthermore, an optimization 
of ammonia flow during the nitridation stage made an 
extraordinary improvement of the InN film's flatness; 
it means that based on the In bilayer controlling depos- 
ition, a moderate, stable, and slow nitridation process 
also plays the key role in growing better-quality InN 
film. Meanwhile, the biaxial strain of InN film was 
gradually relaxing when the parameters of growth 
was optimizing, implying that the mismatch stress of 
InN heteroepitaxy can be well relaxed via this growth 
method. 



Table 1 XRD peak position of (002) diffraction and main 
lattice constants of InN films for our samples 

Sample A Sample B Sample C Sample E Sample F 



lnN(002) (°) 


15.82 


15.83 


15.95 


16.15 


16.19 


c(A) 


5.68 


5.67 


5.63 


5.57 


5.56 


InN(IOI) (°) 


16.65 


16.60 


16.53 


16.43 


16.37 


d101 (A) 


2.70 


2.71 


2.72 


2.73 


2.74 


o(k) 


3.54 


3.56 


3.58 


3.61 


3.63 
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